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An interesting and satisfactory fluid model has been proposed in literature for the the descrip- 
tion of relativistic electron beams. It was obtained with 14 independent variables by imposing the 

O ■ 

entropy principle and the relativity principle. Here the case is considered with an arbitrary number 

*^ ' 

Qh' of independent variables, still satisfying the above mentioned two principles; these lead to conditions 
43 ' whose general solution is here found. We think that the results satisfy also a certain ordering with 



respect to a smallness parameter e measuring the dispersion of the velocity about the mean; this 
ordering generalizes that appearing in literature for the 14 moments case. 



x. 

1 Introduction 

Amendt and Weitzner proposed in [lj and [2] the following system of quasi-linear partial differential 
equations 

d a V a = 0, d a T°P = eF^Va, d a A a ^ = 2eT^F^ a , (1) 

to describe the behavior of relativistic electron beams, where V a is the particle flux vector, T a/3 
is the energy-momentum tensor , F@ a is the electromagnetic field tensor, e is the electron charge 
and A 01 ^ represents the tensor of fluxes. Let's consider the counterpart of the above variables in 
statistical mechanics. They are defined as moments of the distribution function f(x a ,p a ), 

V a = J fp a dP, = J fp a pPdP, A a ^ = J fp a pPp"<dP, (2) 

where p a is the four-momentum of the particle so that we have p a p a = —m 2 and dP = \/—g dpid p^ dp3 
is the invariant element of the momentum space; m is the particle mass. Prom eqs. (|2|) the following 
"trace condition" holds: 

A^g/s-y = -m 2 V a . (3) 
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The closure proposed by Amendt and Weitzner in [T] and [2] is covariant, complete and satisfies the 
above mentioned ordering with respect to e; the resulting system is hyperbolic, but it satisfies the 
trace condition only approximately. In [3] and [3] a new closure has been found. It is very satisfactory 
because it satisfies exactly the condition Q and leads to a symmetric hyperbolic system, but imposes 
the entropy principle only up to second order with respect to thermodynamical equilibrium. In [5] 
the same result has been achieved but up to whatever order. Here the exact general solution for 
the many moments case is found, satisfying all these conditions up to whatever order. Extension 
to very many moments are needed in order to improve on the results of ordinary thermodynamics, 
as shown in [6], page 197. Moreover, we will exploit general equations so that the results may be 
applied also to other materials even if different from electron beams. 

For the case with many moments we have to choose an even number M and an odd number N; after 
that the equations are 

Q J s aa\...a M — Jol\...ol m 

(4) 

d a B aai '" aN = j a i--- a N 

Equations involving lower order tensors are already included in @ because of the following trace 
conditions ([5]). 

All the tensors appearing in the above equations are symmetric and M+N is odd in order to obtain 
independent equations. The counterparts of these variables in statistical mechanics are 

A a ai ...a M = J f p <x p <xi . .. p a M d p^ B a ai ...a N = J yyy*i . . . p <*N d p^ 

from which the trace conditions follows 

2 i«ai...» M _ 2 _ Aaai...a M 2 r>aa 1 ...a N _ 2 _ r>aai...a N /r\ 

— "I Si — H 9aM-lOi M i ~ rn D — D ya N - 1 a N - {<J) 

Let us define the maximum trace of a tensor as the trace of the trace ... of the trace of this tensor, 
so many times as possible. The maximum traces of I ai — a M anc [ Q f J a ^--- a N are ze ro, so that the 
maximum traces of eqs. @ are the conservation laws of mass and of momentum-energy. 
Now there are less independent components in the eqs. (jU) than in the variables i mi, " ftM and 
B aai — aN , so that relations between these variables are needed; aim to find them is called "clo- 
sure problem". We investigate this problem following a commonly used procedure in Extended 
Thermodynamics, i.e., to impose the supplementary conservation law 

d a h a = a > 0, (6) 

that must hold true for all the solutions of the system As a bonus we will obtain a symmetric 
hyperbolic system with nice mathematical properties, such as well-posedness of the Cauchy problem 
and continuous dependence on the initial data (see [7]). By using a property proved by Friedrichs, 
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Lax and Liu [8], [9], [TO], we obtain that eq. © amounts in assuming the existence of parameters 
A Ql ...a M e fj, ai ... aN , called Lagrange Multipliers, such that 

dh a = A ai ... Qm dA aai — aM + Mai ... ajv dB^i-^, (7) 

\ rai...a M , rai...ajv > n 

/ Vii-qm j "i - A^ai '-ajv — u ' 

where /i a is the entropy- (entropy flux) 4- vector. 

We follow now an idea developed by Boillat, Ruggeri and Strumia which they applied in different 
physical contexts concerning both the classical and the relativistic case |12| . |13| . In other 
words, we define 

h '° = - h <* + \ ai ... aM A aa ^ + i^ ai ... aN B a ^- aN , (8) 
and take the Lagrange Multipliers as independent variables. In this way eq. ([7J)i becomes 

dh = A 1 M dX ai ... aM + B 1 N dii ai ... aN , 

from which 

Bh' a 8h' a 
d\ ai -dM d{i ai ... aN 
In this way the tensors appearing in the balance equations (JH) are found as functions of the pa- 
rameters X ai ... aM and i-L ai ...a N i called also mean field, as soon as h' a is known. So it remains to 
find the exact general expression for h a such that both members of eq. ([9]) are symmetric; some 
steps in this direction will be done in section 3. It will be necessary to distinguish between the 
cases N=l and N > 1; the first of these will be analyzed in section 5, the other in section 6. Their 
treatment needs some useful theorems which will be proved in section 4. In the next section we will 
consider the definition of equilibrium and the implication of the above conditions in this state. At 
last conclusions will be drown. 

We consider necessary to mention also the important context of which the present paper becomes a 
part, i.e., Extended Thermodynamics. To this end we conclude this section by citing the pioneering 
papers [13] for the classical context and [15] for the relativistic one, and also the important book [6] 
which contains many results in this context. We reccomand to reed also the important paper [16] 
by Geroch and Lindbloom. 

2 Definition of equilibrium and properties 

Equilibrium is defined as the state described by the independent variables A and \x a such that 

A/3 r ../3 M = A# (/3l/ g 2 • • • 9/3 M -iMi~ m ) 2 ' 
m~0N = V(f3i9foi3 3 ■■■9p N . 1 p N ){-m 2 Y E ^ 1 , (10) 
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from which it follows 



_ n {M 1)!! „aiQ2 . . . „«M-1<*M/ m 2 Y 



'(M + 2)!! r " aMi 
iV!! o iv-i 

{N + 3)!! ^ aai '" a N-i9 9 ' "> i - ■ l- IJ ) 



The physical meaning of this definition is evident when we substitute eqs. (]10p into eq. ©i; by 
using also the trace conditions ([5]), eq. (J7J)i becomes 



d/i a = \dA a + fi ai dB aai 

which is still eq. ([7])i , but in the case M=0, N=l, i.e. in the case we consider only the conservation 
laws of mass and of momentum-energy. 

It is also evident that eqs. (11 1H are identities if M=0, N=l. (Use firstly the identity (M — 1)!! = 

(Af+1)H \ 
M+i >■ 

Eqs. ([9]) now become 



8h a Bh a 



d\ ' d/j, ai 

and the symmetry condition on the second of these is surely satisfied; in fact, from the representation 
theorems [7J, [8] we have that 

h' a = H(\, 1 )v a (12) 



with 7 = yj— fi a fj, a and H (A, 7) is an arbitrary function. It follows 

A " = -8X (13) 

1 f) JT 

B aax = —fi a ^ + Hg aai (14) 

7 07 

which is surely symmetric. 

We stress now that the function H(\,j) has to be arbitrary, if we want that our equations may be 
applied to all materials. In fact, eqs. (fT2"j) - ([Hj) and eq. (JSj) yield 

h a = -nsu a , A a = nu a , B afi = eu a vP + ph afi (15) 



with 



u a = ^, h alS = g a/3 + it V (projector) , 

7 

dH 

n = 1~q^ (particle density), p = H (pressure), 
dH 

e = —H — j—— (energy density), s = — A — J-^jj (entropy density). (16) 

7 ~d\ 
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From (|16p it follows d = ~ds which, compared with the Gibbs relation 

Tds = d (^j + pd (17) 

identifies 7 as with T the absolute temperature. Now (II6P /1 can be used to change variables from 
7, A to 7, p; by substituting its solution A = A(7,p) into eqs. (fT6jh ,5 we obtain the state functions 
n = 71(7, p) and e = e(7,p). Vice versa, if these state functions are assigned, by substituting p from 
(fTB]) ^ into (HS])3,5, these become 

dH 1 r , 

— = -n[ 7 ^(A,7)], 

| r = -_ e [ 7 ,J(A )7 )]--H(A 1 7) (18) 

which are differential equations for the determination of the function H(X, r y). 
Therefore, if we want that our field equations may be used for all materials, i.e., for all possible 
state functions n = n(^,p) and e = e(j,p), then the function H(X,j) must be arbitrary. 
Note also that the integrability condition for eqs. (fT8|) is (e + p)n p — jn^ = ne p . But this is not a 
new condition on the state function because it is the same integrability condition for the equations 

ds d / e \ d I ' 1 



dp ^ dp \n) dp \n 
ds d / e\ ^ ( 1 

dj ^ dj\n) dj \n / 

which are equivalent to the Gibbs relation reported above in eq. (1171) . 

We conclude this section showing that an expression for h' a at equilibrium is 



h a = J F(\,p u p»)p a dp (19) 
where the function F(X,Y) is related to the distribution function at equilibrium by 

-^F(XX) = fe q (X,Y). (20) 



The expression ([19]) will be useful in the sequel and is equivalent to (|T2|) with 

H(\,j) = m 3 / F(X, jm cosh p) sink 2 p dp. (21) 

7 Jo 

which gives the relation between H and F. 

It follows that F(X,Y) is arbitrary, because i?(A, 7) is arbitrary, and (|19p is the most general 

expression for h a at equilibrium. A particular case follows from (|20p and from 

w 

feq(X,Y) = x+y 

e~ ± 1 

which is the Juttner |19j . |20j distribution function at equilibrium, where k is the Boltzmann con- 
stant, the upper and lower signs refer to Fermions and Bosons, respectively, h is the Plank's constant 
and w is equal to 2s+l for particles with spin 1^. 
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3 The entropy principle 

To impose eqs. ([9]) we have to find the most general expression of h' a such that the left hand sides 
are symmetric. We will refer to this as "the symmetry condition". An exact particular solution of 
this condition is 

h[ a = [ F(X ai ... aM p a \..p a ^,^ 1 ...^K..p^)p a dP, 



as it can be easily verified, where the function F(X,Y) has been determined in the previous section 
(eqs. (|18p — (|2ip ) in terms of the state functions at equilibrium. This solution is more general than 
the corresponding one in the kinetic approach |21j . where the particular case F(X, Y) = F(X + Y) 
is considered, but it is not still the most general one. We aim here to find this most general solution. 
To this end let us note that the symmetry condition for eqs. ([9]) reads 

dh'^ a dh'l a 

0, T\ = 0; 



d^a.\\a.2—OiM ^Mc*i]«2"-c*iv 

if we subtract from these their expressions with hf instead of h' a , we find 



0, = 0. (22) 



with Ah' a = h' a - h\ a . 

But in the previous section we have found that {ti{*) eq is the most general expression for (h' a ) eq , so 
that we have 

(Ah' a ) eq = 0. (23) 

In this way, we have now to find the most general solution of eqs. (|22p and (|23p . after that we will 
have 

h' a = h[ a + Ah' a . (24) 
To exploit eqs. (|22|) and (|23p . let's firstly consider the Taylor expansion of Ah' a around equilibrium: 

00 i 

Ah ' a = E m^'"^'"^---^^---^ (25) 

h,k=0 

where the multiindex notation Ai = ■ ■ ■ cti M , Bj = aj 1 ■ ■ ■ aj N has been used, and, moreover 

, .n.i| ••• i;, /; • •/>', _ [ d h+k Ah' a \ 

h ' k \dX Al —dX Ah dfAB 1 —d f JLB k ) { ' 

\ / eq. 

and 

' ~ i\i 
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denote the deviation of the Lagrange multipliers from their value (llOh at equilibrium. Note also 
that their traces are zero, as consequence of eq. (fTT|) , 

Because of the symmetry shown by eqs. ()22[) it is possible to exchange the index a with each other 
index taken from those included in A{ or Bj and moreover each A{ can be exchanged with each 
other A s or B r . So the tensor A'fju is symmetric with respect to every couple of indexes. 

Let's consider 9Xa ...qx^q^ ■■■<9mb. d e P encun S on ^71— 7m as a composite function through A 7l ... 7M 
and A; after that let us take its derivative with respect to A 7l ... 7M and calculate the result at 
equilibrium. We obtain 

p. r ,otA\"-A h B 1 "-B h ~ p. ri aA\--A h B\--B k 
r aAv-A h+1 B v -B k _ OU h,k °*Pi~Pm , OU h,k OX _ 



p )( ^aA 1 ---A h B 1 ---B k . 



(71 . . . „7m) _ (7172 . . . a lM~HM) aFl „ . . . „ fl 2^— — I + 



p. ( ~,aA 1 -A h B 1 -B k 

<3±Kk (7172 . . 7A/ . 17M)2 M-l !! 2 f 

9A y y (M + l)!! v 7 ' v ; 



where (fTTjb and (|2T|) i have been used. 
Multiplying both sides by <? 7l72 ■ • • 9^ M -hm we obtain 

F)r° 

#7172 yiM-iiM — k m ) ■ 

Similarly, let's consider ^ — ® x + ^ — T77 q- — depending on ^ 771 ... 7JV1 as a composite function trough 
M771-7N-1 an d A*/35 after that let us take its derivative with respect to /U 771 ... 7JV _ 1 and calculate the 
result at equilibrium. We obtain 

p ir iaA i —A h B i —B k „~ n r< aA 1 -A h B 1 -B k _ 
c ,aAi-A fc Bi-B fc+1 = ^fe.fc OUPPi-Pn-! OU h,k C^ff = 

-Pn-i / i 77i'"7iV-i 

/i 771 ... 7jv _ 1 

- ^ fl (7 (I 71 • • • fl 7 "- 1 ' - (7 7172 • ■ ■ fl 7 "- 21 "- 1 ^ «•••()« a ^8 — 1 + 

~~ a^...^ V ^ ^ Av ~ 1 (/3 2 Pn - 2Pn ~ i] {n + J 

ariOiAi-A h Bi-B k 

+ ZZhh ,(7 q 7172 o7j y-27iV- 1 ) 8 Nl] ( m 2)^ 

+ d H 9 P 9 9 8 (/Y + 3)!! ( m > > 

where i and (|27|) -? have been used. 

Multiplying both sides by <7 7l72 • • • 5 , 7JV _ 2 7iv-i we obtain 

aA 1 ...A h B 1 ...B fc /3 7r .. 7jv _ 1 _ ° U h,k (-m 2 )^ (30] 

U /l,fe+l #7l72 g~iN-2~iN-l — Q^^ \ 171 ) \ 6U ) 



So, from eq. (|26j) we have obtained the compatibility conditions (j29j) and (|30|) . 

Let's proof the vice versa, i.e., that eq. (|26|) is a consequence of all the other equations. To this end 
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it is firstly useful to show a property of Ah a . If we take its derivative with respect to A 7l ... 7Af we 
obtain, with passages like that used in eq. ([28]) . that 



dAh' a ~ 1 dC%f- A ^- B * dX 

h,k=0 n ,M 

Eh r ,aAx-A h B 1 -B k 'X 7 ~ ~ (71 7m) ^ 

Tm c h,k k *A 1 ---\A h ^B 1 ---VB k gg 1 1 ---gZf > 



where we have used eq. (|29jh So we have proved that derivation of eq. (|25|) with respect to A 7l ... 7M 
is equivalent to its derivation with respect to A 7l ... 7M . 

With analogous passages, but by using eq. ([30]) . it is possible to prove that derivation of eq. (f25j) 
with respect to /i 7l ... 7JV is equivalent to its derivation with respect to /I 7l ... 7jv . 

We are now ready to prove that eq. ([26]) is a consequence of the other equations. To this end, let 
us take the h th derivative of eq. (|25p with respect to A 7l ... 7M , then its k th derivative with respect to 
|U 71 ... W and calculate the result at equilibrium; by using eqs. (f29|) . ([30]) and the above mentioned 
property , we obtain eq. ([26]) . In other words we can forget eq. ([26]) and retain only eqs. ([29]) and 

m. 

So it remains to solve eqs. (j29p and (|30f) in the unknown symmetric tensors C' h ' k . 
We notice that both sides in equation (|29[) and the left-hand side of eq. (|30p are symmetric; then 
we have to impose that the right hand side of this last one is also symmetric. In other words, both 
the tensors C' h ' k and their derivatives with respect to up are symmetric. The tensors satisfying this 
property are elements of a family T which will be characterized in the following section. Moreover, 
interesting properties of this family T will be shown and they are useful to exploit our eqs. (|29p 
and (|30p . The effective exploitation will be accomplished in section 5 for the case N = 1 and in 
section 6 for the case N > 1. We complete the present section simply reporting the results, so that 
they are available to whom is not interested in the proofs. 
For the case N > 1, they are 



C^ 1 k ' aMh+kN+1 = (jh,kg(aia 2 . . . ga 2a -ia 2s ^ 2s +i . . . ^a Mh+Nk+1 ) ^ ^-g 



s=0 
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with 



„[ Mh+Nk+l l Tfcl ,„ f Mh+Nk+1 1 | 

C fe > fe = 2 2 • 5 " 



s! (M/i + iVfc + l-2s)! 

Mh+k(N-l)~2 

- e -M h ~(N + i )k+ 2s V- [Mfe + fc(JV-l) + l + 2[|]]!! 

4^ [hM + k(N - 1) - 2q - 2]\\ 



q=0 

Mh+(N+l)k 



and c k>q = c k>q (X) restricted only by c k ^ q = c fc _ lj9 for q = 0, • • • , [ Mfe +( fc - 1 K Ar - 1 )- 2 l . 
For the case iV = 1, and consequently k=0, the results are 

Mh 
2 

(J^ l ~" a Mh+l _ (jhg(aia2 . . . g02s-ia2s ^a-23+1 . . . ^a M h+i) ^33^ 



with 



s=0 



Mfa-2 

r h _ 2 Mfe-2s ( i l^) ! 1 -6-Afh+2 S V- (Mfe + 1)!! 

s! (M/t + l-2s)! 7 4^ (Mh-2q-2)\\ 



and ^ « = c fe q(A) restricted only by g = c^-i g for g = 0, • 



[M(/t-l)-2] 



4 On the family T and its properties 

Let us characterize now the family T of tensorial functions of a scalar A and of a time-like 4- vector 
Hp, which are symmetric together with their derivative with respect to ug. 
We will prove now that 

Proposition 1 Each element of the family T can be written as 

[f] 

s=0 

urai/i i/ie scalars ^satisfying the condition 



1 



57 7 



2s + (n-2s + 2)(n-2s + l)0™_ 1 = /or s = 1, • • • - . (36) 



n 



(Note that among the terms 0™ we may call leading term the one with the highest value of s, i.e. 
4>J n y Once the leading term is known, we can find all the other terms present in eq. (1351) thanks to 
eg.' ([2|j. 
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Proof. From the representation theorems we know that eq. (|35h is the most general expression of 
a symmetric tensorial function depending on A and /ig. Taking into account that = — the 
derivative of (p ai '" a ™(A, jip) with respect to jip is 



Lai— an 



dfip 



E 



h n ,fi 



$7 7 



M (aia2 n a2s-ia2s ,,«2s+i ,, Q »0 _i_ 

y ' ■ y p ' ' ' Li t 



(37) 



To be symmetric, the expression above must be equal to its symmetric part with respect to 
«i • • • a n (3, i.e. to 



,^2,-1^^+1...^^) + 



d0" 1 ( ai q 2 

s=0 ' ' 

[§]+! 

+ ^ 0"_ iy (aiQ2 • • • ^-l«2S^2S + l . . . ^«n/)( n -25 + 2) 



5=1 



$7 7 



gOl a 2 . . . _|_ 



+ E 



c?7 7 



+ (n - 2s + 2)C_i 



1 



g(«l a 2 . . . _|_ 



+ 



E 



«? 1 



re + 1 - 2s 
n + 1 



c^7 7 



- + (n - 2s + 2)#_! 



2s 
n + 1 



gP( a l . . . g a 2s-20t2s-l ^Ot2s . . . ^«„) _|_ 



g( a l a 2 . . . ^02 S -102 S ^02 S +1 . . . ^%)^ 



where in the second row we have substituted s = S — 1. In the third row we have reported the term 
coming from the first row with s=0, and that from the second row with S = [§ ] + 1; in the fourth 
row the remaining terms from first and second row. 
Comparing the result with eq. (|37p . we obtain 



E_0<£_s_*_ ( ai a 2 n a 2s -ia2s ,,oi2s+i „««) 
9 77 9 9 ^ M 



s=0 



c?7 7 



+ E 



s=l 



(?7 7 



+ (re -2s + 2) 



n + 1 - 2s 
n + 1 



„(ai"2 „Q!2s-l02s ,,a2s + l 

y y a* Z- 4 
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and 



J2$(n - 2s)gP {ai ■ ■ ■ g a ^»2s+i . . . 



s=0 

<^ n 1 n - 1 
[2 1 ' 



n 
L2 



gOl"2 . . . gQn)/3 _|_ 



[»]-> 

+ £ 



s=0 



$7 7 



+ (n-2s)C 



2s + 2 
n + 1 



where we have considered s = S + 1 and then S = s in the last term. 
So we have 



1 



and 



d-y 7 



2s + (n - 2s + 2)(n - 2s + l)$?_i = for s = 1, • 
2s + 2) + (n - 2s) (n - 2s - = for s = 0, • • 



n - 2 



$7 7 

We can see that the second of these equations coincides with the first one. So the characteristic 
condition for T is the above reported eq. ([36]) . 
This completes the proof of Proposition [TJ ■ 

It will be useful for the sequel to note that, thanks to (|36j) . eq. (|37l) becomes 



C +1 (A,7)5 (Q1Q2 ... 5 Q 2 S -ia 2s/i a 2s +l 



(38) 



with 



1 ggffi 

7 87 



±1" 



(39) 



/£+i = «±i( n _2s + 2)«^_ 1 fors = l,---,[^ 
This allows to prove the following 

Proposition 2 J/ g it follows that also ^ G J 7 . 

Proof. Because of eq. ([38]) , to demonstrate the theorem it will be sufficient to prove that eq. ([36 
holds true also with n+1 instead of n, i.e. 

1_ , . . _ . . Tn + r 



S7 7 

For s=l, thanks to (|39p . it becomes 



2s + (n - 2s + 3)(n - 2s + 2)<#JT 1 1 = for s = 1 



(40) 



37 7 2 



fn + 1) + (ra+ l)n 
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that is identically satisfied. 

Instead, for s = 2, • • • , [^^1 , thanks to ([39|b . it becomes 

2s n + 1 ^ _ ^ + 2) + ^ _ ^ + ^ _ ^ + ^ _ 2g + = Q 



that is eq. (|55|) with s-1 instead of s, after having divided it for 2s-2. ■ 
It is also important the following 

Proposition 3 I/^ 1 "^ 1 £ T, then 6 a ^- a ^ £ T exists such that (h a ^- a «+^ = 2±2 . 

Proof. In fact the integrability condition for the problem 

,ai-a n+ i _ U V 



^+1 

is 

= 



which is surely satisfied because £ J? and, consequently, it and its derivative with respect 

to fJ- Qn+1 are symmetric. 

We note also that, if (p ai "' an is a particular solution of this problem, then the general one is 



tai— an 



0ai-o„ + 0(A) ff («i Q 2 . . . g a„-ia n ) j f n ig 

if n is odd, 



with (j)(X) a scalar function. 



Proposition 4 If 6 ai " a "-+ r £ T, then <A"i- tt » £ T exists such that 6 a i~ a »+r -- 

Moreover, if <f) ai "' an is a particular solution of this problem, then the general one is 

; ai ...a n + ^4 V] ^(A) 9 (°i«2 • • • g a "+^ a « +2i ^ an+l ■ ■ ■ fi an+2i if n is even, 



Lai ■■■a r , 



Jai-an + ^| = 2 J ^(A) 5 («i^ • • • g a ^ a " +2i+1 ^a n+ i ■ ■ ■ Ma„ +2i+1 n *s odd, 
wft 0i(A) scalar functions. 

Proof. We can prove this proposition with the iterative procedure. 
It holds for r=l for the Proposition [3l 

Let us assume now that it holds for r = r and prove that it is satisfied also when r = r + 1. 

If ^ ) Q i'" Q! "+~+i £ J? we can apply this proposition with n+1 instead of n and r instead of r. Then we 

have that d> cll '" oln + 1 £ T exists such that 
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But, for the first part of Proposition[3]we have that Ap- 1 '"^ £ ^"exists such that (h a i--- Q n+i -_ 

°h L a n + 1 

which, substituted in (JH]) gives 



^£>l-«n+r+l __ 

^Mon+i " " " dfla n+ - +1 

So the existence of solutions has been proved. 

If Ql '" a " is a particular of these solutions, we have 



(42) 



which, together with eq. ([HZ]) implies 



d 



iai— a n itii-a, 



0. 



By applying the second part of this Proposition 01 we conclude that 







if n+1 is even, 



(«1«2 . . . n «n + 2i + l"n + 2i+2) ,, . . . ,, 

i/ rOn+2 ^Q n+ 2i + 2 



EU J 0i(A)5 
if n+1 is odd. 

This relation can be integrated and gives 

r ~-i i _ 

2 ' (b (\)—^— n( Q i a 2 . . . n a n +2ia n +2i+l) ,. ..... 

if n is odd, 

[ r-2 1 _ 

2 ' (b (\)—^— tf( a l a 2 . . . „«„+2 l +l«n+2i+2) ,. ..... i 

+</>o(A)^ aia2 • • • if n is even. 

So we have obtained the second part of this proposition, with f + 1 instead of r, and 



&(A) 



6 i _ 1 (A)^ if n is even (For this case we have put i=I-l) 



This completes the proof of Proposition [U ■ 

Proposition 5 // </, Q i" a ™ e T, then 0«i" a ™+2 g JF exists such that 



^- a " +2 y an+ian+2 = r n 



Moreover, its leading term is 



h n+2 _,,-2(n+3-[^l) + 2) 



^n^ 3 -^^ + /^(A) 



(43) 



(44) 



wit/i (/>"„-| leading term of (f> ai "' an and an arbitrary function. 

2~ 
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Proof. Because of (j) a ^'" a n+2 g j: ^ ^ j^g ^he f orm (j35|) with coefficients satisfying eq. (|36p and n+2 
instead of n, i.e., 



mi— a„+2 , 



! (A,^) = J] C +2 (A, 7 )5 {QlQ2 ---5 Q2s - ia2 ^ a2a+1 ---^" +2) 

with 



m+2 1 

-2s + (n-2s + 4)(n-2s + 3)^+ 1 2 = for s = 1, 



n + 2 



s 7 7 - • • • < 45 > 

Let's substitute this expression of 4> ai ' " an + 2 in the left hand side of eq. (|43p and explicit the 
symmetrization, so obtaining 

m 

f n+2 1 



E ^" +2 ( A ' 7) , ^T^rW " 1 " 2 • • • 5 a2s - 3 ° 2s - 2 ^ 2s - 1 • • • ^ + 
^ (n + 2){n + 1) 



a=l 

fn+21 



[n + 2){n + 1) 



3=1 



^ 2 ^ (n + 2-2 S )(n + 2-2 S -l) 
hi (A ' 7) (n + 2)(n + l) 

Blending the first two sums and putting S=s-1 we obtain 



.2 

iQi-a„+2 



V (.^+2, (^ + 1)^-^ + 2) 2 2 (n + 2-2s)(n + l-2s) ] 

g« n+1 * n+ 2 - Z>\^+i 4 („ + 2)(n + l) 1 7j (n + 2)(n + l) J 

g{aia2 . . . ^a2s-ia2s^ i a2s+i . . . 

With the use of this expression and of eq. ([33]) for ai '" Qn , eq. ([35]) becomes 

M f ( 8 + i)( w - 8 + 2) +2 _ 2 (n + 2-2s)(n + l-2s) ) 

^l^ +l4 (n + 2)(n+l) + (7) (n + 2)(n + l) /' 

[f] 

g("i"2 . . . ^Q2 s -ia2 s ^a2s+i . . . ^ a n) — ^rig{oL\cx2 , . . g02s-i02s^02 s +i . . . ^o n) 

i.e. 

€ = ( n + 2)(n+l) ^ 4{S + " S + 2) + ^ +2 (- 7 ')( n + 2 - 2 S )(n + 1 - 2,)] , (46) 
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fora = 0,-.- ,[§]. 

If we use eqs. ([36]) and (14"51) . this expression becomes 



J s+l 



1 



^7 (n + 2)(n + l) 



in+2 d6 n+2 

s+2 -4(s + 2)(n - 8 + 2) + — 7 2 )(n - 2s)(n - 1 - 2s) 



(47) 



that is eq. (|46p with s+1 instead of s, derivated with respect to 7 and with another use of eq. (|45j) . 
So we have proved that if eq. (|46p holds true for a particular value of s, it will hold true also for all 
lower values of s; so it suffices to impose eq. (j4*6j) for the bigger value of s, i.e., 

^ 2 [f]+ 2 



r [fJ (n + 2)(n + l) 

where eq. ([431) has been used. The general solution of this equation is reported in eq. 
It will be useful also the following 



(48) 



then the leading term of (f> ai "' an g an _ 



in—2r 

5 rn-2rl 



(2 [s±i] -2r-l)!! 

(2m-i)n 



/(A)7 



-2(3+p), 



is ^ 








2 


5a n -lo 




"n + r 






2 



2r - 2 - 2p, 2 



n+ 1 



4 - 2p) (49) 



where, if a < h, rj(a, b) denotes the product of all even numbers between a and b, while if a = b + r 
then rj(a, b) = 1. 



Proof. Let us prove eq. (|49j) with the iterative procedure. 
In the case r=0 it is an identity. 

Let us assume that eq. (|49|) holds up to the integer r and let us prove it with r+1 instead of r. 
Let us distinguish the cases with n odd and with n even; for the first one, by applying eq. (|4"8j) with 
n-2r-2 instead of n, we obtain 



,ra-2r-2 

■> ' n-1r-Z 



1 



n — 2r 



in -LV + :i)CL +1^-<K~IU 
— 2 — 07 — 2 — 



(n-2r-2)!! 
nJl 



f(X)r](n — 1 — 2p — 2r, n — 3 — 2p)7 



-2(3+p), 



n - 2r + 3 - 6 - 2p) 



that is eq. (1491) with r+1 instead of r. 

In the case with n even, by applying eq. (|48p with n-2r-2 instead of n, we obtain 



m-2r-2 

- ) n-2r-2 



1 



n — 2r — 1 



(n - 2r + 2) 



in— It 1 „ 
9„-2r + 7 



C?7 



in—2r 

' n — 2r 



(n-2r-3)!! 



(n- 1)!! 

that is eq. (I49p with r+1 instead of r. 



f(X)n(n - 2 - 2p - 2r, n - 4 - 2p)7~ 2(3+p) (n - 2r - 4 - 2p) 
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Proposition 7 If Q i'" Q »™ g jr with leading term /(A)7 2 ( 3+p ), with p a non negative integer 
such that p < m — \Jj\ — 1 or p > m — + r — 2, i/ien a i -- Q m+2r g eX 2 S i s gyc/j that 

Aai---a m a m+ ia m+ 2---a m+ 2r-ia m +2r n . . . n — A,ai—a m 

Y !Ja m + ia m+ 2 i/Om + 2r-lO!m + 2r V 

Moreover the leading term of (p ai '" arn+2r is 

(m + 2r)! (2[f])H (2m - 2 [g ] - 2 P - 4) !! 2{3+p) , f w - 2(s+B ^- m ) 

m! (2 [f ] + 2r) !! (2m - 2 [f ] + 2r - 2p - 4) !! /(A)7 + 1. -M^ 

(50) 

u>ii/i fi r (X) arbitrary functions. 

Proof. Let us prove this proposition with the iterative procedure. 
It is easy to verify that eq. (|50p holds for r=0. 

Let us assume that it holds up to the integer r and prove that it holds also with r+1 instead of r. 
Then we have to face the problem 

jai-a m a m+ ia m+ 2-«m+2r+iam+2r+2 „ „ Ap-i---a m 

V ya m+1 a m+2 ya m+2r +ia m +2r+2 — V \ OL ) 

By defining 

J.ai---Ctm+2r iai-a m +2 r am + 2r + l«m+2r+2„ ( 

9 — 9 ya m + 2r + ia m+ 2r+2> 

the problem (|5ip becomes that of the present proposition, which we have assumed holding. There- 
fore, the leading term of Q i'" a ™.+2>- i s f|5Q[) and it remains to face only the problem f)52[) : by applying 
Proposition [5] with n=m+2r we find that the leading term of <^ Q i" a m+2r+2 ^ s 

" (m + 2r)\ (2[f])H 



2 (3+m+2r-[lg±^]) ( m + 2r + l )( m + ^ + 2) 



2 r m+2?-+L' i 



ml 



(2 [f]+2r)! 



r-1 

[ro+2r+21 



(2m 2 [™] 2p 4)!! 2 p+2m+4r-2| g±|^ ]-l , \ ^ f . m -24+2^1,/-, , fm+2r+2 /n 

(2m-2[f]+2r-2 P -4)!! /(A)7 + ^ /v( " *T + 



7 



i=0 

(m + 2r + l)(m + 2r + 2) (m + 2r)! (2 [f ])!! (2m - 2 [f ] - 2p - 4)!! 

2 [ m+2r+2 j m! (2 [f] +2r)!! (2m- 2 [s] +2r-2p-4)! 



, m 1 ^-2 P ,STf m (m + 2r + l)(m + 2r + 2) 

• / ^_2p + 2m + 4r-2[^^± 2 f± 2 ] 7 ' 2 [™±f±2] 

1 „-2i-2m-6+ar^ai , «n+2r+2 / ^ n (m + 2r + l)(m + 2r + 2) _ 2 f m+2r+3 _[m±^±2i') 

^7 + /^m+2r+2 j V-^J ^±J±I] ^ 

that is eq. (|50p with r+1 instead of r and 



i,r+l 



, (m+2r+l)(m+2r+2) 1 , . n 

Ji,r 2 ^ m +2r+2 j 1r-1i ' ' ' 

f m+2r+2 m (m+2r+l)(m+2r+2) , . 
J r m+gr+2 j l A J 2 JHi+2r+2j — 



This completes the proof. ■ 

We conclude this section with the 
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Proposition 8 If n is an even number, the tensor g( aia2 ■ ■ ■ g an ~ ian ' fj, an _ r+1 ■ ■ ■ fj, a „ belongs to T 
and, moreover, 



n — r 
[ 2 J 



n (aiot2 . . . „a n _ia„) 

1/ !J H'Oln-r+l H>a r , 



lU—T (aiQ!2 . . . n a2s-ia2s . ."2S + 1 . . . ,«n-r) 



^ ?r,..' 9 ' " - ■■■(, - • -'ii -■■ • ■ ■// 



s=0 



(53) 



m-r 



, r!(n-r)! (_>v 2 V 

(2s+2r-ra)!!(2s)!!(n-r-2s)!(ra-l)!! V ' / 







/or r < 1 and s = \J L 2 L ~\ , 

/or r < 1 and < s < - 1, 

/or r > 2 and § -r < s < [^], 

/or r > 2 and < s < § — r — 1. 



(54) 



Proof. It is easy to note that the above tensor is an element of T because it and its derivative with 

respect to (ip are manifestly symmetric. 

Let us prove (|54|) with the iterative procedure. 

It is easy to verify that it holds when r=0. 

When r=l it is a consequence of 



„("1«2 rfiln-lCI-n) . . _ „(Ql«2 „«i.-l)«i.„ _„("lQf2 „ttn-3«n-2„«n-l) 

y y H'an — y y H-a n — y y p 



(55) 



When r=2, by using eq. (|55j) we have 



>l a 2... rt a »-l Q »)„ 

^(«1"2 . . . ^an-3«n-2)^ttn-l _|_ ( n _ 2)g a "- 1 ( Ql . . . ^ a n-4«n-3 ^Qn-2 ) 

n — 1 L 

_~,2 „ o 

_J n {aia2 . . . n a n -3«n-2) _|_ 



-.9 



•3 



5 aia2 ---ff an - SQn -V an - 3 Pa n _ 2 ) 



n — V ' n — 1' 

from which (|53p and (|54p with r=2. 

Let us assume now that our proposition holds up to a fixed integer r and let us prove that it holds 
also with r+1 instead of r. By multiplying eq. (|53p times fi a „_ r we find 



n ( a l a 2 . . . A-lOn) ,, ..... 
y y P«n-r H'a-n 



^slhr^-^y^ 2 ■ ■ ■ g a2S ~ ia2S fi a2S+1 ■ ■ ■ M a «— + 

s=o 



n — r 



n — r 
I 2 J 



s=0 



H r 2£^ (Q1Q2 _ _ _ a2s _ ia2 a2s+1 _ _ _ an-r-j)/ 2x 

n — r 
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where, in the first summation, we have put s=S+l. 
It follows (1531) with r+1 instead of r and 



Lfl—1 1 

9 s,r+l 



n—r 2s+2 .2 in— r n—r— 2s 



s+l,r n—r 



7 



,r n—r 



for s = 0) ... ,[«=!] -1, 

and for s = [^y^] in the case r even, 



I ts,t n—r 



for s 



n—i — 1 



in the case r odd. 



We have taken into account that n-r-2s=0 when s = and r is even. 

By using this result and eq. (|5^|) we have that 



• For < s < § — r — 2, both </>™ + [ and 4>™ r are zero, from which ^"r+i 1 



0. 



in-i — l 



(n-r-l)! 



For s = § - r - 1, we have 4> n s / = and y, s>r+1 - (n _ 2r _ 2 )!!(n-i)!! • 
For s = 2 — r,"- , ~ ^ anc ^ ^ or s = iT~\ m ^ ne case r even we have 



in— r— 1 



(r + 1)! (n-r-l)! 



(2s + 2r + 2 - n)!! (2s)!! (n-r-l - 2s)! (n - 1)!! 



2\s+r+l- 



For s = - 2 1 and r odd we have (j^r+i 1 



(r)! (n-r-l)! 



(r-1)!! (n-r-l)!! (n-1)!! ( 7 ) 2 • 

In this way we have obtained (|54j) with r+1 instead of r as we desired to prove. ■ 



5 The case N=l. 

When N=l only equation (|29|) with k=0 has to be exploited. 

We prove now that it amounts to giving the following expression for the leading term of C aAl "' Ah 

Mh-2 

M h d h c hjq (X) ( 1 V (Mh + l)\\ 



9=0 ' \ * / 



with c h , q = Cfc-i^ for q = 0, • 



i 1 ) (Mh-2q-2) 

M(h-l) -2 
2 



and C(l „ fa , = M , . . , are arbitrary functions of A. (50) 

Let us prove this with the iterative procedure. 

It holds for h=0, because in its right hand side there are no terms, and its left hand side is zero for 
eqs. (J23]) and (EH). 
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Let us apply now the Proposition [7] to eq. (129p with k=0, r = -g-, m=Mh+l. We find that 



2 

A/h-2 



r h+i -6 V- r m 2^( fe+ i) ^ +1 ^ g (A) /lV (Mfe + 1)!! [M(h + !) + !]! 

f(h+i) 7 ^ 1 ' dX^ 1 V7 2 / (M/i-2g-2)!! (M/i + 1)! 



M-2 
2 



(Mfe)!! (Mfr-2g-2)H ^ ^-2(3 + Mi+,) 



[M(h + 1)]!! [M(/i + 1) - 2q - 2]! 
which is eq. ([56]) i with h+1 instead of h (c^+i^ = Ch A has to be used), i = q 



Mh 



c h+l,q 



Mh-21 



Ch, q for q = 0, • • • , [ 

[M(fe+l)-2 g -2]!! , 2 \-f(h+l)f* f nr „_Mh M(ft+l)-2 

[M(/i+l)+l]M v m / •> g—MA M. 101 « — 2 2 



with f* mh M defined by 4^if*Mh m_ ~ fn—MA — • 

<f 2 > 2 " y 2 > 2 ^ 2 ' 2 

We note that from eq. it follows 



M(h - j) - 2 

Ch,q = c h -j, q for g = 0, . (57) 



Also this relation can be proved with the iterative procedure. It holds when j=0. From (|57p and 
([55|) 9 it follows c/i j(? = c-h-j,q = c h-j-i,q for g = 0, • • • , M ^~2~^~ 2 that is eq. (|5T|) with j+1 instead 
of j, as we desired to prove. 

Let us search now the other coefficients C^. To this end, it is better to prove firstly that from eq. 
(1361) with n=Mh+l it follows 



h Mh+i = ( ,y g! (Mh + l-2s)\ 8 r Mh+X 
K ' (s-r)\{Mh + l-2s + 2r)\d( n 2 y Vs 



Let us prove this with the iterative procedure. 

It holds for r=0. Let us assume that it holds up to the index r. From eq. (|36p with n=Mh+l and 
s-r instead of s, we find 

iMh+i ~ 2« + 2r 1 d ,Mh+i _ 



J s-r-l 



7 {Mh + 3 - 2s + 2r) {Mh + 2-2s + 2r) <9 7 s_r 
4(s-r) ' 9 ^ Mh+l 



{Mh + 3 - 2s + 2r){Mh + 2 - 2s + 2r) d{-f 2 ) s ~ r 

(-4) 



. r+1 s! {Mh + l-2s)\ ^+1 > Ih+1 



(s-r - 1)! {Mh + 3 - 2s + 2r)\ d{-f 2 ) r+1 Ys 
where in the last passage eq. (|58p has been used. In this way we have proved that (|58|) holds when 
r+l replaces r; this completes the proof of eq. (|58p and we use now it to find Cg. 
If we write eq. ([58]) with s = ^-h and r = s — s*, jointly to eq. (|56|) . we find that 

Mh-2 



c\ = ( _ 4) f^(f ^) ! i y ( _ m 2 } f 

(s*)! {Mh + l-2s*)\ f^ Q { ' 



d h c hs {\) {Mh + 1)\\ , M h _ s *{q + 2 + M-h-s*)\ ( iy + ' i+ - h - 



d\ h {Mh-2-2q)\V ' (9 + 2)! \j 



„2 
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that is, the above mentioned eq. (l34"l) . 



6 The case N > 1 

In this case we have to impose eqs. (|29p and (|30p . But it will be firstly useful to transform them 
into more easy equations. To this end we may use the notation Bi = ■ ■ ■ aj iV „ 1 (similar to the 
already used multindex notation Bi = a,i x ■ ■ ■ a.% m ) and Proposition 2J we obtain that the tensor 
D'^ Al ' AhBl " Bk exists, such that 



r ,uA\---A h fi\B\--fi k B k _ u r .a/i 1 --A h o 1 --D k 

-d Hx -d Hk B ^ (59) 



After that eq. (f30|) becomes 

d k+1 ri aA 1 -A h B 1 -B k71 - 7N _ 1 ^ _ ^_^_^^^_ F) aAi-A h B 1 -B k( _ to 2n^± 

For the Proposition |4] with n=hM+k(N-l)+l the general solution of this equation is 

n aAi-A h B 1 -B kJ1 - 1N ^ 1 _ n aAi-A fc Bi-B fc /_ 2s^i , 

^h,k+l i/7172 " ' ' yiN-2"lN-l — ^h,k V J 

+ £ ^AkW9 [aM --- AhBl --- Bkah ^ (60) 

i=0 

Let us transform now the unknown tensors D' hk into the tensors E' h ' k , according to the following 
rule 

■ ) aAi--A h _B 1 - i? fe _ v aA 1 -A h B 1 -B k 



u h,k - ^h,k + 



i=0 

with ifii t h,k(ty defined with the iterative method by 

ipi,h,a = 0, 



4>i,h,k+l 



[Mh + k(N -l) + 2i + 4]!! [hM + (it + l)(N - 1) + 2i + 1]!! 
[KM + (fc + l)(N - 1) + 2z + 4]!! [hM + k{N - 1) + 2i + 

o iV-1 

^i,h,k{-'m ) 2 +4> ithtk 

k — 2 — 3 

in the case k even and i = 0, • • • , and in the case k odd and i = 0, • • • , 

2 2 

[Af^ + Jfe(JV r -l) + 2i + 4]H [hM+(k + l)(N -l) + 2i + l]\\ 
^' h ' k+1 [hM + (fc + 1)(N - 1) + 2i + 4]!! [/iM + A;(iV - 1) + 2i + 1]!! Afc 

— 1 

in the case k odd and i = . 
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We note that also E' h ' k G T . 

Thanks to this transformation, eq. (|60p becomes 



^jfc+l 57172 • • • 9in-21n-i - h h,k K~ m ) ■ i bZ ) 

Let us note how many 4- vectors \x a intervene in the second term in the right hand side of (|6ip ; they 
are 2i + 1 < 2 [ k -^] + 1 < k — 2 + 1 < k; therefore eq. (|60p . substituted in (159p transforms it into 

r ,aA 1 ---A h (3 1 B 1 ---(3 k B k _ 9 ^aA 1 -A h B 1 -B k 

Eqs. ([62]) and ([63]) substitute now eqs. ([59]) and ([60]) ; consequently, these last one can now be left 
out! 

Eq. ([63]) with k=0 now yields 

C h,0 ~ h h,G ( b4 i 

which is restricted by eq. (|29p with k=0. After that eq. (|62p will give, with an iterative procedure, 
E' h ' k+1 . Let us firstly deduce, from this procedure, some properties. 

Property 1 The leading term of E^ Al "' Ah is multiplied by a polynomial in the variable 4j. 

Proof. This property is evident from eq. 



Property 2 The leading term f E? Al "' AhBl " Bk is multiplied by a polynomial in the variable 



J h,k 

l 

r" 

Proof. Let us prove this with the iterative procedure. 
It holds when k=0, for property 1. 

Let us assume it up to the index k. From eq. (|62p . by applying Proposition [7] with m=Mh+(N- 
l)k+l, r = ^y 1 -, we find that property 2 holds also when k+1 replaces k. 
This completes its proof. ■ 

By substituting now eq. ([63]) into eq. ([29]) . this last one becomes 
dfifr ■ ■ ■ dij,3 k 

For Proposition [5] with n=Mh+(N-l)k+l, r=k, it follows that 



r^aAi-AhBi-Bk-yx^-JM-iJU n _ ( ^2\4f " F aAi-A fe Bi-B fc 

^h+ljk i/7i72 " " " y7M-i7M I d\ ' k 



T?o<A 1 --A h B 1 --B k ~f 1 ~j2--lM-ilM n _ (_ rr ,2\^- 9 w aA 1 -A h B 1 -B k 

^h+hk i/7172 ' ' ' y-lM-HM \ rn J "oT-^i 



h+l,k y7i72 yjM-l7M \ ) Q^h,k ~ 

= J2 <MA)S Ml '--^ Bl '-- B * Q1 '-- a24+1 Va 1 • • • /Wn" (65) 

2=0 

where the notation g a ^ - a ^« = g( a i a 2 . . . g^in-xam) h as Deen used. 

But, for the above property 2, jointly with the Propositions [6] and [8] we have that the left hand 
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side of eq. (|65p has a leading term of the type multiplied by a polynomial in the variable 

while its right hand side has a leading term polynomial in 7 2 . It follows necessarily that both these 



leading terms are zero; in other words, both sides of eq. (|65p are zero. In particular, from the left 
hand side we obtain 

„aAi-A h Bi-B k j l72 -'jM-l'yM n _ i m 2\^- 9 a A 1 -A h B 1 -B k / fifi >. 

^h+l,k #7172 ' 9-yM-lJM - \~ m ) 2 'Q\ tj h,k ■ i bb i 

Therefore, we have to impose only eqs. 

HS2) and dSSD; after that, the tensor c^ Ahlil ^ k defined 
by eq. ([63]) is the more general solution of eq. ([29]) and ([30]) . 

6.1 A consequence of eqs. ( 1621 ) and ( 1661 ) 



We will see now that, as a consequence of eqs. ([62]) and ([66]) . we may find the leading term of 
E^ 1 ' AhBl ' Bk except for a set of arbitrary functions of the single variable A. This result is the 
subsequent ([70]) . 

To this end let us firstly prove that 



(iV-l)fc-2 



t (-^)^V(A) [ ^:^;, [(^- 1 ).-^(1)°. (67) 
where the last N-l functions Ck, q are arbitrary functions of A and the remainder are 

Ck,q = C k -l,q- (68) 

Let us prove it with the iterative procedure. 

It holds when k=0, because in this case the right hand side has no terms, while the left hand side 
is zero for eqs. ([Ml) with h=0, ([23]) and (p5]) . 

^ — N 

2 ' 



Let us apply Proposition [7] to eq. ([62]) with h=0, r = ^-i-, m=k(N-l)+l. We find that 



(JV-l)fc-2 

r,o,k+i - 6 ( 1 V [(iV-l)fc + l]!! UAT 

^ + D = 7 g Cfc 4^J [(N-l)k-2 q ]!! [{N - 1)k - 2q] - 

, _2^(fc+l) [(N-l)k + l + N-l}\ [(N - l)k]l\ [(N - l)k -2q- 2]!! 

1 ' [(JV - l)jfe + 1]! [(JV-l)jfe + JV-l]U[(JV-l)(A; + l)-2g-2]n 

V-3 
2 

-[6+k(N-l)+2i] 



j=0 

Then we have found eq. (|67[) with k+1 instead of k, i = q — k^^- and 



c k+l,q 



Ck, q for q = 0, • • • , (JV 2 )fc 2 , 



Jg-kE^l^ [(JV-l)(A:+l)+lj! J 2 ' tor g 



[(AT-l)(fc+l)-2g-2]!! ( _ 2 \-^(k+l) fm . „ _ (N-l)k (jV-l)(fc+l)-2 

2 ' ' 2 
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This completes the proof. 

It will be useful in the sequel also the 

Property 3 c ks = c k -i, g for q = 0, • • ■ C^zlMzllzg , 

Proof. Let us prove this with the iterative procedure. 

It is obvious when i=0. 

From eq. (|68l) with k-i instead of k, we have 

^k,q Ck—i,q Ck—i — l,q 

for q = 0, • • • — 1 ^ fc 2 1 ^ 2 . So the property is valid also when i+1 replaces i. 



Now, from eq. (|62p it follows 

T?aAx-A h Bx-B k _ p fiiAi-A h B 1 -B k B k+1 B k+2 , 2 \-(iV-l) 

^ft,* -^,fc+2 9 Bk+1 9 Bk+2 { m ) 

and so on, until 

-^h,k+Mh 9 Bk+1 ---9 Bk+Mh { m ) ■ 



By applying h times eq. (|66j) . the above equation becomes 

w aA v -A h B v -B k _ d h aA 1 ...A h B 1 ...B k7v .^ Mh{N _ 2) {N-2)Mh 
^h^k ~ Qyh 0,k+Mh i/7i72 " ' i/7Mh(JV-2)-i7ii/h(JV-2) V " L > ■ 

(69) 

In this way all the E'^\ are determined in terms of E' ' k . 

From eqs. ([59]). (|67|) and Proposition [6] with n=Mh+(N-l)k+l+Mh(N-2) and r = 4f h(iV - 2) we 
find that 



(JV-l)(fc+Mfa)~2 
2 



p h,k _ -6 \^ / ^^k+^^k+Mh^ / 1 \ 

[M/t + (iV- !)& + !]!! 



^ifc+Mh^Cfe+M^g /IV [(JV-l)(fc + M/») + l]l! 



d\ h Vt 2 / [(iV-l)(fc + M^)-2g-2]I! 



[(M/i + fc)(iV -!) + !]!! 



?/[Jfe(JV - 1) + M/i - 2q, (Mh + k)(N - 1) - 2 - 2g]. 



This expression can be simplified because, for g > Mh+k ( N 1 ) we have M/i + [N — l)k — 2q < 0, 



(M/i + k)(N — 1) — 2 — 2q > 0, so that at least one, among the factors intervening in ry(- • • ,•••), is 
zero. Then we can restrict to the values with q < Mh + k (N-i)-2 anc i ) consequently, 

(N-l)k+Mh-2 

w h,k _ _ 6 , 2^ H « d h c fc+JfM / 1 y [Mh + (N-l)k + 1]!! 

Mh+( f- 1)fc 7 ^ 1 j dA^ V7V [(JV-l)fc + Mfe-2g-2]H" 

(70) 

This agrees with the above used property 2. 
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6.2 Equivalence of eqs. ([62]) and JBBj) with eq. ( TTDT) 

The result ([70]) has been proved as a consequence of eqs. ([62]) and ([66]) , We prove now the vice 
versa, i.e., that eqs. (|52"j) and ([55]) become identities when eq. ([TO"]) is used. 

Let us begin with eq. (|62p : For the Proposition [6] with n=hM+(k+l)(N-l)+l and r = ^r-, and 
([70]) , the leading term of the left hand side of eq. ([62]) is 

(Af-l)(fc+l) + Mh-2 
2 



-6 v [(iV-l)(fc + l)+Mfe + l]!! 2 N^± (k+1)+ M* d h c k+1+ Mh, q 

P5 [(iV-l)(fc + l)+M/i-2g-2]!! 1 j dA^ 

[(N — l)k + Mh + l]\ 



7V [(JV-l)(fc+l)+M/i+l]!l 



?](M/i + fe(JV - 1) - 2g, Af/» + (fc + 1)(N - 1) - 2 - 2g). 



If g > Mfc+fc 2 (7V ~ 1) we have KM + fc(JV - 1) - 2q < 0, hM + (k + 1)(N - 1) - 2 - 2q > 0, so 
that r/(- • • , • • • ) = 0; therefore we can limit to values with a < Mh + k ( N ~ 1 ) _ x and the above 
leading term becomes equal to the right hand side of eq. ([70]) premultiplied by (— m 2 )^-^- (use the 
property ([68]) with k+l+Mh instead of k, i.e., Ck+i+Mh,q = Ck+Mh,q for g = 0, ■ • • , Mfe+fc (^ f ~ 1 )~ 2 < 
{N - l){k + I + Mh) - 2). 
Therefore, eq. ([62]) is an identity. 

Let us now prove the same thing for (|66[) : For the Proposition \5\ with n=M(h+l)+k(N-l)+l and 
r = 4^, and ([70]) . the leading term of the left hand side of eq. ([66]) is 

(JV-l)fc+M(fe+l)-2 

_ 6 [(N- l)k + M(h + 1) + 1]!! 2 jv^fc+ MWj ^ +1 c fc+Mfe+M , g 

7 ^ [(N -l)k + M(h + l)-2q-2]\r m) 2 dX^ 1 

1 \ 1 \(]\T _ l A/f/, _|_ 

[ (^ + aW)+V (M '' + - » - 2 '' M <* + D + *(^-D-a- 2,). 

When q > Mh+k ( N 1 ) we h aV e ry(- • • , • • • ) = so that there remain terms with q < Mh+k (^ l) 2 ^ 
which is the right hand side of eq. ([70]) multiplied by (— m 2 )~2~ and derived with respect to A (use 
Ck+Mh+M,q = Ck+Mh,q for (j = 0, ■ • • , Mh+k(N l) 2 h y s f or property 3 written with i=M 

and k+Mh+M instead of k because Mh + k(N - 1) - 2 < (N - l){k + Mh) - 2). Therefore, eq. 
([66]) is an identity. 



6.3 Determination of the tensor c%f'" AflBv " Bk 



After having imposed eqs. (]62p and (|66|) . let us now impose eq. (j63p for the determination of the 
tensor (j^v A hBi-B k rp Q ^ s en( \^ us firstly note that from eq. (f39]) 9 with n odd and s = 
we find that the leading term of ^ is 

^1 = ^1=1- (71) 

~2~ 2 
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From eq. (j39H ? with n+1 instead of n and s = we find that the leading term of a ^ is 

(h n+2_ l? n + 2 +1 n + 2 8 
~~5~ n + 1 7 c*7 ""s - 

where in the last passage eq. (I71D and eq. (I36p with n+1 instead of n and s = have been used. 
If <fin-i depends on 7 by means of 7 2 , it follows that also depends on 7 by means of ■j 2 and is 

^ = -2^ + 2)-^^ (72) 

It follows that 

(In fact, it holds for r=0. Let us assume that it holds up to an integer r. Eq. (|72p with n+2r instead 
of n becomes 

_ ^ 2(n + 2r + 2) A_^-2 :i = (_ 2) r+l (« + ^ + 2 ) !! ^ + ' , 



n _ 2, i — j — . (73) 



where ([73]) has been used. The result is again ([73j) but with r+1 instead of r; this completes the 
proof). 

From eq. (|7ip . with n+2r instead of n, and for eq. (|72|) . it follows 

m+2r+l _ / \r (" + 2r ) !! in 

This result and eq. (f72|) give 



(-2)[l] 



; ,(,i + 2 [§])!! alt] 



l-a- J n!! 9( 7 2)[aJ 2 



This relation, with n=Mh+k(N-l)+l, and eq. (|70p allows to obtain from eq. (|63p that the leading 



term of (j'*M---A h B-L---B k j g 



Mfe+fc(iV-l)-2 
2 



d h c k+MKq {q + 2+ [§])! / 1 ^ 



dA h (? + 2)! V7 2 

This result allows to determine the other coefficients C^' fc . 

To this end, it will be useful to note firstly that, from eq. (|36p with n=Mh+Nk+l it follows 



(74) 



s\ (Mh + Nk + l-2s)\ d r 
(s - rj\ {Mh + Nk + 1 - 2s + 2r)\ <9( 7 2 
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(In fact, this relation holds for r=0. Let us assume that it also holds up to an index r; from eq. 
(|36p with n=Mh+Nk+l and s-r instead of s we find that 

r = -2s + 2r 1 d_ = 

Ys ~ r - 1 7 (Mh + N k + 3 - 2s + 2r) (Mh + Nk + 2 - 2s + 2r) <9 7 s ~ r 

1 d ,„ 



^ (M/t + iVfc + 3 - 2s + 2r)(Mh + Nk + 2- 2s + 2r) <9( 7 2 ) 
( 4)r+ l * ! (Mfr + iV fc + l-2s)! g^ 1 xn 



(s-r- 1)! (M/t + Nk + 3 - 2s + 2r)\ d(^ 2 ) r+1 s 

where in the last passage eq. (j75|) has been used. The result is eq. ()75|) . with r+1 instead of r, and 
this completes its proof.). 

Eq. (J75J) with s = [ Mfe+ ^ fc+1 ] and r = s - s*, jointly with eq. (J7U) shows that 

M = [MMp±l] | ( Af/l + Nk + i _ 2 [Mfej^hi] ), 

1 J s*! (M/t + iVfc + l-2s*)! 

mh + k(N~-l)-2 

[|] A (Mfe + fc(JY-l) + l-2[§])! , y-u , m 
W ^ (M/i + fc(iV - 1) -2g-2)! 1 ' ' 



d h c k+Mh , q (<i+2+ m ( n r^i-- (g+ 2 + [§] + ggg] -**)! 

dA" ( g + 2)! 1 ) ' (q + 2 + [§])! 



i v+3+[!M 



i_ | • Uh+JVfc+l j 



from which the above reported eq. f)32[) . taking into account that 

"M/i + iV/c + l 



Mh + Nk + 1-2 



0! = 1 if k is odd, 
l! = 1 if k is even, 



and moreover, that [Mh±Nk±lj + [|] = Mfe +(f+ 1 ) fc , 



7 Conclusions. 

We consider the present results very satisfactory, because the closure which is found constitutes 
a consistent improvement of the previous one appearing in literature. It takes into account an 
arbitrary, but fixed, number of moments; moreover it satisfies the constraints exactly and satisfies 
also the supplementary conservation law related to the entropy principle. This fact guarantees the 
existence and uniqueness, well-posedness and stability of the solution of the initial value problem. 
Moreover it allows also further investigations, such as that concerning the warm plasma ordering 
and that concerning the behavior of the subsystems. These may be the subject of future works. 
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